Structure of the core of magnetic vortices in d-wave superconductors with a 
subdominant triplet pairing mechanism 
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The quasiparticle states found in the vortex core of a high-T c cuprate superconductor may be 
probed by scanning tunneling spectroscopy. Results of such experiments have revealed typical 
spectra that are quite different from what is seen in conventional low-Tc superconductors. In 
particular the Caroli-deGennes-Matricon state at E ~ in the core center is not seen. Instead, in a 
high-Tc vortex core, quasiparticle states are found at energies that are at a sizable fraction of the gap 
energy. One explanation for this could be that a finite amplitude of a competing order parameter 
stabilizes in the vortex-core center. Here I will explore the possibility of nucleating a vortex-core 
state that locally breaks inversion symmetry. The vortex-core order parameter is of mixed parity, 
a [d + zp]-wave, and the quasiparticle spectra in the core center lacks the E — states. 

PACS numbers: 74.25.Ha, 74.50.+r, 74.55,+v 



I. INTRODUCTION 

In conventional, nominally clean, type-II superconduc- 
tors the quasiparticle spectrum in a vortex core was de- 
scribed by Caroli, deGennes, and MatriconP They found 
that quasiparticle states are localized in the vortex core, 
and that these same states carry the currents that screen 
the magnetic flux line from penetrating the interior of 
the superconductor. A direct measurement of vortex- 
core states was done by scanning tunneling spectroscopy 
(STS) 2 and the measured data could be quantitatively 
explained by theory. 3 Performing STS measurements on 
vortices in high-T c cuprates, on YBCO 4 or on BSCCO 5 7 
(see also references in Ref. [8]) revealed a very differ- 
ent generic picture; the vortex core in a high-T c super- 
conductor does not harbor pronounced core states. By 
theory^MI these core states should be present also for an 
order parameter of d-wave symmetry relevant for high- 
Tc superconductors and thus readily be seen in STS. In- 
stead, non-dispersing vortex states are seen at energies 
corresponding to ~ 30% of bulk value of the supercon- 
ducting energy gap, Ao. 

Theoretical suggestions to explain the vortex-core 
spectra seen in high-T c cuprates include possi- 
ble subdominant singlet-paring amplitudes JHHUD anti- 
ferromagnetic order stabilizing in the vortex corepEU 
the normal state pseudogap phase made visible in the 
coreP^, or competition between d-density wave and d- 
wave superconducting order In this paper I will explore 
the possibility of nucleating a vortex-core state of mixed 
parity and show that this state may have a finite ampli- 
tude of a p-wave order parameter in the vortex center. 
The possibility of a mixed singlet-triplet pairing state in 
a vortex core has been suggested to occur in an s-wave 
superconductor as a result of spin-orbit coupling.^ In 
the context of a d-wave superconductors both spin-orbit 
coupling and the presence of a Zeeman coupling has been 
considered, finding a singlet-triplet mixing in the vortex 



phase in the Ginzburg-Landau region (0 « T < j^fEHHl 
Here, I will use the microscopic quasiclassical theory to 
show, given an attractive pairing interaction in a p-wave 
channel and a weak intrinsic Zeeman coupling to the 
magnetic field to break the spin-singlet symmetry of the 
parent d-wave supeconductor, that a sizable p-wave or- 
der parameter may stabilize locally in the d-wave vor- 
tex core region. Using this self- consistently determined 
order-parameter field I then compute the spatially re- 
solved local density of states in the vortex. It turns out 
that the quasiparticle spectra seen in the STS on the 
high-T c cuprates^sj can to large extent be reproduced 
theoretically as a direct consequence of this triplet su- 
perconducting core order. 

The scenario is that a finite amplitude of a compet- 
ing order parameter stabilizes in the vortex-core. An 
isolated singly-quantized vortex has the asymptotic or- 
der parameter A (12) = A^e 1 ^ as one circles the core. 
The phase winding of 2tt corresponds to a center-of-mass 
angular momentum, L c z m A(R) = f -^A(R) = HA(R), 
of the Cooper pairs. Approaching the vortex core, the 
order-parameter amplitude is gradually suppressed and 
vanishes in the origin of the vortex so to maintain single- 
valuedness of A (12). To have a finite order-parameter 
amplitude in the core the center-of-mass angular momen- 
tum can rotate in to an internal orbital angular momen- 
tum of the Cooper pair, L° z rh A(p,R) = |^|~A(p, 12) = 

hA(p, 12). 22 This scenario occurs in the B-phase of su- 

perfluid 3 He where A-phase and double-core states are 

found to be energetically favorable to a normal-state core 

in different regions of the pressure-temperature phase 
diagramPEUM] 

A d-wave superconductor has singlet-pairing symme- 
try and to have a finite order parameter in the core with 
L° rb = h a p-wave order parameter is needed in the 
vortex-core region. This requires; i) an attractive pairing 
interaction in a triplet channel, ii) a symmetry breaking 
field that introduce a seed of a triplet component. The 
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FIG. 1: [Color online] A schematic picture of the vortex-core state of a d-wave superconductor in the presence of a pairing 
attraction in a triplet channel is shown in panel a). The core state [R < £o) may have a superconducting core with a sizable 
order-parameter amplitude with a p x ± 2p y -symmetry. In panel b) the coupling constants computed using a spin-fluctuation 
mediated pairing interaction equation M and a simplified generic band structure of the cuprates. The parameter 8 tunes from 
predominantly anti-ferromagnetic (S — 0) to ferromagnetic (S = tt) spin- fluctuations. The width of the curves signifies the span 
from an overdoped (|Ar(<5 = 0)| smaller) to an underdoped (|Ar(£ = 0)| bigger) material. Attraction in the spin-triplet pairing 
channel is sizable for all 5. 



first condition I argue to be intrinsic in a spin-fluctuation 
mediated pairing relevant for high-T c superconductors 
(see Ref. 27 and references therein). This interaction 
can support condensation into both spin-singlet and spin- 
triplet superconducting states!^ The second condition is 
readily given by the weak Zeeman coupling to an exter- 
nal field always present in a vortex. This is in particular 
true for extreme type-II superconductors where penetra- 
tion depth is far larger that the coherence length. In 
this paper I neglect orbital effects and the screening of 
the external magnetic field and assume that the external 
magnetic field is constant over the vortex-core region. 



II. MODEL 

For highly anisotropic spin-fluctuations, \ z = Xzz ^ 
Xxx,yy = X ± i a susceptibility that can be tuned from pre- 
dominantly antiferromagnetic (5 « 0) to ferromagnetic 
(6 ~ 7r) spin-fluctuations can be modeled as 
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Xq is the overall amplitude and £ s fi is the spin-spin corre- 
lation length which is typically a few lattice constants (a) 
in the cuprates. Introducing a simple coupling g between 
the spin fluctuations and the quasiparticles the pairing 
interaction is V(p — p') = V(q) = Mfg 2 x z {o) — X z {q)i 
Aff being the total density of states at the Fermi level. 
The resulting gap equation allows for three channels of 
pairing, one spin-singlet with V s (q) = X z (q)/^ an d two 



spin-triplet channels, one with V tz (q) = X z (q)/2 having 
the d- vector parallel to z and one Vt±(q) = — x 2 (g)/2 for 
which d _L z. Using equation ([!]) the gap equation in a 
weak-coupling approximation reads 

AM = -T E {VM~li)n{ti)fM\Zn)) n . (2) 

\e n \<e c 

Here f x (p^; e n ), (x = s,tz,tl.), is the anomalous prop- 
agator at Matsubara frequency e n and momentum p, 
and nip,) = Mi^TV^^)! -1 )]^- 1 use a lm " 
earized version of eq. assuming the factorization 

Vxfe-li) = Er ^rG^r^), together with a tight- 
binding parameterization of the band structure relevant 
for BSCCOp^to compute the eigenvalue spectra for the 
possible pairing symmetries as a function of 5 keeping 
£ 3 fi = 2a. To each eigenvalue Ap belongs a set of ba- 
sis functions ^rfe) which may be classified according to 
the irreducible representations (T) of the crystal group 
D^. The resulting eigenvalues as function of doping of 
BSCCO and the degree of incommensuration of the spin 
fluctuations are shown in figure [l]b. Attractive eigenval- 
ues (Ar < 0) are found for the even-parity representa- 
tions Aig^g, B\g,2g with strongest attraction in the B\ Q - 
channel (d x i _ y i -wave) followed by the ^^-channel (ex- 
tended s-wave) for dominantly anti-ferromagnetic spin 
fluctuations. Also the odd-parity representation E u has 
attractive eigenvalues. Here it is the channel with d \\ z 
that is most attractive, closest in value to that of the 
I^-channel. This parallel to what was found earlier in 
the case of Sr 2 Ru04pHl 

Next, I need a theory to self-consistently solve for a 
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FIG. 2: [Color online] The order-parameter amplitudes computed at T = 0.05T C with Xe u = 0.9\b 19 and h = 0.02Ao, i.e. 
2% of the zero-temperature d-wave gap, Ao. The asymptotic d-wave orderparameter having a phase winding of 2tt is as seen 
suppressed in the core and heals to its bulk amplitude over roughly 5 £o away from the core. In the vortex core a substantial 
triplet orderparameter ~ p x -i y is stabilized filling an area of ~ 7r£o with an amplitude close to 50% of the bulk value of the 
d-wave gap. This component has no phase winding but relative phase of —tt/2 to the d-wave component. Away from the core 
a second component ~ p x +i y with a phase winding of 4ty appears. 



vortex structure in the presence of a Zeeman field and 
competing order-parameter symmetries and use the qua- 
siclassical theory, a leading order theory in A/Ep(<^. 1), 
as presented in e.g. Ref. [22 j. The central object is the 
Green's function g(j^,R;e) — (1/a) f d^ p rsG(p 1 R, e), 
obeying the Eilenberger equation 

^V£(ft, R; e) + [s%-h&-A(%,R),g{&,R- e)} = 0, (3) 



and normalization condition g(ft,i2;6:) 2 



The 



quasiclassical " Hamiltonian" , e%— h • a — A (ft, R), is 
a 4 x 4 matrix in combined particle-hole (7^=1,2,3) and 
spin space {(Ji= x ,y,z)- The order-parameter matrix in (|3| 
reads 



A(a,fl) 



A(ft,H) 
A(ft,fl) 



(4) 



where A(ft,fl) = [A s (ft,#) + A* (ft, R) • cr]ia y is a 
2x2 spin-matrix order parameter. The superscripts 
refer to spin-singlet (s) and spin-triplet (t) compo- 
nents of the order parameter. Particle-hole components 
are related via the "tilde" -symmetry a(fy,R]£,t) = 
%,R] —£*,t) with * denoting complex conjugation. 
This gives A(ft,fl) = ia y [A s *(-&, R) - A** (—ft, R) • 
a] = ia y [A s *(fa,R) + A (j^,R)-a]. The Zeeman term 



in Q reads h • & = diag[/i • tr, — a y h • cra y ] with the Zee- 
man field pinning the spin-quantization axis to the z-axis, 
h er = —\ibBg z . h is assumed to be homogeneous and 
small, \h\ = h <C |Ao|, to weakly break parity (Ao is the 
zero-temperature d-wave gap). 

The Eilenberger equation Q may be solved by intro- 
ducing the following two spin-matrix coherence functions 
7 = (7s+7t -cr)ia y and 7 = ia y ( ; y s —'j t -a) parametrizing 
the Retarded quasiclassical Green's function} 10 * 22 1 30 1 



-inN 



1 + 77 
-27 
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-I-77 



9 f 
f 9 



(5) 



with N = diag[(l — 77) _1 , (1 — 77) -1 ]. The Advanced 
function is given as g A = T$g R ^T% and the Matsubara 
function as g M (e n ) = g R (e + iO ie n ). The mixed- 
parity orderparameter components are linear combina- 
tions of the singlet and the z-component of the triplet 
part as A±(ft,fl) = ±[A s (ft,#) ± A£(ft,fl)] (and 
A± (ft , R) = ± [ A s * (ft , R) ± A'/ (ft , H)] ) . This leads, to- 
gether with the Zeeman-shifted frequency e± = e±[1bB, 
to a separation in to two pseudo-spin bands (±) with 
different orderparameters A±(ft,i2) and equation Q is 



4 



written as two pairs of scalar Riccati equations 



iv F -\/j± + 2£ ± 7 ± 
ivp- V7± — 2e±j± 



- 7± A ± 7 ± 
-7 ± A ± 7 ± 



(6) 
(7) 



one for each spin-band. The two equations are solved 
by numerical integration along straight lines, or trajec- 
tories, s(x) = Sq for 7± and j± as described 
in e.g. Ref. [22 . Once 7±(^, R] e) and 7±(^, -R; e) are 
obtained, the order-parameter fields are calculated using 
^ as 



a5(h) = —\B lg T (ykM) 

\e n \<e c 



(/+"/-) 



^^(h) = —\e u t J2 (yk,±m {J± i L1 )0 

\e n \<e c 

where /± = f±(p F1 R]e n ) are the anomalous functions 
on spin-band ±. {' ' ')p F = J ^ s tne avera S e over the 
momentum direction on the Fermi surface, with the 
angle <p p giving the angle the momentum makes to 
the crystal a-axis (x-axis in the figures). In the vortex 
calculations the conventional basis functions yB lg (ib) = 
\/2cos20p and yE u ,±{%) = -\/2(cos0p ± zsin0 p ) are 
used. The paring interaction and cut-off frequency, e c , 
are eliminated in favor of the transition temperature T c as 
-\-\ g = mT/T c + En>o( n + V 2 ) -1 - The subdominant 
interaction Xe u is introduced in the self-consistent calcu- 
lations as fraction of the dominant one and is treated as 
a parameter free to explore. 



III. RESULTS 



Equations ( 6ff ) and ( 8|9 ) are iterated until self- 
consistency is reached. In figure [2j the structure of 
a mixed-parity vortex is shown. A substantial triplet 
p-wave order parameter may be nucleated in the sin- 
glet d-wave vortex core with both p-wave components 
p x ± ip y present. The p-wave component, A p , with 

r orb 

z,p—wave 

and it carries no phase winding. The 

A ■ tttHVi T orb 



^Td-wave wm finite m the core center 



p-wave compo- 
nent, A p /, with L° r °,_ wave = ~L c z m d _ wave has a finite 
amplitude on the phase boundary separating the singlet 
and triplet order parameters, at R ~ £o = hvp/^T c1 
around which its phase winds by 4tt. This so that 

L Zp'-wave + LC Z'-wave = L d-wave' A11 amplitudes re- 
tain the four- fold symmetry of the d x 2_ y 2 amplitude as 
seen in the contour plots displayed in figure [2] 

The nucleation of a p-wave order parameter is depen- 
dent on a finite Zeeman field. In figure [3] the ampli- 
tude A p (0) is displayed as a function of h for different 
strength of Xe u - A p (0) is finite for all coupling strengths 
and grows with increasing Zeeman field. For larger cou- 
plings, Xe u > 0.7A# lg , the onset of A p (0) at small fields 
becomes increasingly nonlinear with sharp onset of the 
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FIG. 3: [Color online] The amplitude of the p-wave order pa- 
rameter in the vortex core, A p (0), computed at T — 0.05T C , is 
shown in panel a) as a function of the symmetry-breaking field 
h for various ratios Xe u /Xb 19 - For large ratios, Xe u /Xb 19 < 
0.9, the p-wave stabilizes even as h — » (but h finite). The 
temperature dependence of A p (0) is shown in panel b) at 
h = O.OOlAo and at h = 0.02 A . For the larger field A p (0; T) 
is finite at higher temperatures and grows with decreasing 
temperature while for h = O.OOlAo, A p (0;T) has a distinct 
temperature below which it grows rapidly to its low-T value. 
For smaller ratios, Xe u /Xb 19 < 0.7, the p-wave core phase is 
close to linearly dependent on h. 



sub-dominant order parameter at h > 0. In panel b) of 
figure [3] the temperature dependence of A p (0) is shown. 
While the p-wave amplitude is finite at all temperatures 
there is a transition from a field induced triplet order pa- 
rameter at larger fields, when h > 0.02Ao, and at high 
temperatures T > 0.1T C , to an intrinsic phase transition 
at vanishingly small fields h < O.OIAq and low temper- 
atures. For weaker triplet pairing strength the intrinsic 
phase transition in to a d + ip vortex state is pushed to 
lower temperatures. 

The emergence of a mixed-parity state in the vor- 
tex center at low temperatures has a profound ef- 
fect on the quasiparticle spectra. In general the spa- 
tially resolved density of states (DoS), NT t(R',e) = 
— ^Im(Tr[fsg R (iir, R',e + 20 + )])^, will show evidence of 
the sub-gap Caroli-deGennes-Matricon states that carry 
the screening current of a vortex.^^ Resolving the spec- 
tra also in position on the Fermi surface one finds 
for a pure d-wave vortex, on trajectories tangential to 
the asymptotic phase winding, that the bound states 
in the core have the qualitative dispersion with distance 
or 'impact parameter' b from the vortex-core E(j^,b) ~ 
±| A(jj,)| tanh(6/6o)- Here, bo ~ £0 is a scale factor and 
± denotes if the momentum direction is parallel (+) or 
anti-parallel (-) to the phase winding of the asymptotic 
order parameter A(j^)e^ at b <C 0. Introducing an 
imaginary p-wave order parameter in the core splits the 
quasiparticle spectra in to two branches, one for each 
spin band a(= ±). This is displayed in figure [4] The 
qualitative quasiparticle dispersion of the core states is 
now modified as E a (f^^b) « ±|A(j^)| tanh[(6 =p s a b)/bo] 
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FIG. 4: [Color online] The trajectory-resolved DoS 
N±(jfr,R m j€) evaluated on the cut along the y-axis through 
the origin. The spectra are taken in the indicated points sep- 
arated by 0.5^o (see the inserts). The phase winding of the 
d-wave is counter clock wise in this case. The direction of 
% is chosen (anti-)parallel to the x-axis in (lower) upper pair 
of panels. The presence of a subdominant order-parameter 
amplitude introduces a shift of the the state in the core cen- 
ter to a finite energy (this state is marked by an arrow in 
each panel). This energy shift is negative (positive) for the 
branch N+(N-). There still exist zero-energy states on both 
branches, but these are now found away from the core center. 

where s± = ±1. The offset, 6, is a direct consequence 
of a finite order parameter in the vortex core center and 
shifts the zero-energy state to a finite impact parame- 
ter away from the core center. In the core center the 
quasiparticle state is shifted to finite energy, £^(^,0) ~ 
— s a \A(fy)\ tanh(6/&o)- Note that both states, with and 
against the phase winding (±) on one branch are shifted 
to the same energy. This leads to a suppression of the 
screening current in the core area, \R\ < l£o- 

The features of the trajectory resolved DoS are de- 
tectable in the trajectory averaged total DoS, N Tot (R; e), 
which is directly related to the tunneling conductance 
measured by STSP In figure [H] the DoS calculated at 
T = 0.025T C with \ Eu = 0.9X Blg and h = 0.02 A is 
shown. N Tot (R;e) on a ray through the vortex core 
along an anti-node (panel a) lacks a zero energy state 
(ZES) in the core center. The core state is pushed to 
E a (0) ~ ^0.3Ao and this state have very little disper- 
sion with position on the ray. ZES are found at a dis- 
tance b ~ ±0.5^o from the core center with half the spec- 
tral weight of the ZES in a pure d-wave vortex as the 



two spin bands are shifted differently by A p (0). The 
features in the DoS are also generally broadened by the 
angle average (• • • )p^. The corresponding spin-band re- 
solved DoS is shown in panel c for spin band (-). On a 
ray through the vortex core along a node (panel b) the 
DoS also lack ZES in the core center but the core state 
^a(0) ~ t0.3A has more of a dispersion with small b 
compared to that in the anti-nodal direction. This is due 
to the linear opening of the energy gap around the node, 
I Ad (0)1 ~ \(dA/d(j))(6 - (j) n ode)\ of the d-wave gap. In 
panels (e-i) in figure [5] the spatially resolved DoS is dis- 
played on a 6£o x 6£o-square with center in the vortex core 
at different fixed energies. The ZES form ring around 
the core center and the cores states, E a (0) w ^0.3Ao, 
extends along the anti- nodes. The overall shape of the 
vortex is a doubling of the star-shaped DoS found in the 
pure d-wave vortex coreP^ The doubling shows up as a 
square lattice in the DoS amplitude, and the lattice con- 
stant is set by the magnitude of the induced triplet order 
parameter, A p (0). 



IV. CONCLUIONS 

In closing, I have shown that a mixed-parity d + ip- 
vortex state is possible to stabilize in a high-T c super- 
conductor. The weak Zeeman coupling to the exter- 
nal magnetic field gives a sufficient seed to nucleate a 
p-wave order parameter in the vortex core. The nec- 
essary attractive triplet-pairing channel is supported by 
spin-fluctuation mediated pairing, argued relevant for the 
cuprates. This new core state is directly detectable in 
STS measurements of the quasiparticle spectra and I find, 
within the limits of quasiclassical theory, good agreement 
with existing experimental dataP"^ A further study using 
the Bogoliubov-de Gennes equations is needed to quan- 
titatively compare the theory with experiments. This 
as the quasiclassical theory does not self-consistently re- 
solve the angular-momentum quantization of the Caroli- 
deGennes-Matricon states. This quantization gives a fi- 
nite shift, or a mini gap, of the lowest energy state from 
the Fermi surface ~ A 2 /i?pP^For the high-T c cuprates 
this mini gap may be sizable as A/Ep ^0.1 and may 
very well remove the states found at E = 0.0 away from 
the core center. 
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discussing results, has been extremely valuable. 



1 C. Caroli, P. G. deGennes, and J. Matricon, Phys. Lett. 9, 3 F. Gygi and M. Schliiter, Phys. Rev. B. 43, 7609 (1991) 
307 (1964) 4 I. Maggio-Aprile, C. Renner, A. Erb, E. Walker, and 

2 H. F. Hess, R. B Robinson, R. Dynes, J. Valles Jr and 0. Fischer, Phys. Rev. Lett. 75, 2754 (1995) 

J. V. Waszczak, Phys. Rev. Lett. 62, 214 (1989) 5 S. H. Pan, E. W. Hudson, A. K. Gupta, K. W. Ng, 



6 
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N-(R) e) for one spin-band. There is a dispersion of the core states as a function of distance from the core but as seen in the 
spin-resolved DoS, the total DoS also consists of two non-dispersing states with large spectral weight located at e = ±0.2Ao. 
In panels (e-i) a spatial map of the DoS is shown at various e. 
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